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The characteristic predictions of chiral sohton models - the Skyrme model and its extentions - 
, are discussed. The chiral soliton models predictions of low-lying dibaryon states qualitatively 

agree with recent evidence for the existence of narrow dibaryons in reactions of inelastic proton 
scattering on deuterons and double photon radiation pp ppjj- The connection between magnetic 
' momentum operators and tensors of inertia valid for arbitrary SU{2) skyrmion configurations allows 

to estimate the electromagnetic decay width of some states of interest. Another kind of predictions 
are multibaryons with nontrivial flavour - strangeness, charm or bottom, which can be found, in 
, particular, in high energy heavy ions collisions. It is shown that the large B multiskyrmions given 

' by rational map ansaetze can be described within the domain wall approximation, or as spherical 

bubble with energy and baryon number density concentrated at its boundary. 

^ ■ 

m ■ 1 Introduction 
cn 

2 I The chiral sohton approach provides a very economical method of description of baryonic 
systems with different baryon numbers, starting with few basic concepts and ingredients 
incorporated in the model lagrangian [jl|, The latter is truncated lagrangian of effective 
I field theories widely used for description and explanation of low-energy meson and baryon 
interactions 0]. Baryons or baryonic systems appear within this approach as quantized 
solitonic solutions of equations of motion, characterized by the so called winding number 
^ , or topological charge. If the concept of topological soliton models is accepted and baryons 
are skyrmions indeed, then it is clear why there is isospin in the Nature: the number of 
generators of the SU{2) isospin group, 3, coincides with the number of the space dimensions, 
thus allowing for the correlation between SU{2) chiral fields and space coordinates resulting 
in appearence of topological solitons. 

As it was found numerically, the chiral field configurations of lowest energy possess 
different topological properties - the shape of the mass and B-number distribution - for 
different values of B. It is a sphere for B = 1 hedgehog , a torus for i? = 2 |^ , tetrahedron 
for B = 3, cube for B = 4 0, and higher polyhedrons for greater baryon numbers [0, §|, [^l- 
A paradoxical feature of the whole approach is that the baryons/nucleons individuality is 
absent in the lowest energy static configurations (note, that any of known lowest energy 
configurations can be made of a number of slightly deformed tori). It is believed that 
the standard picture of nuclei should appear when non-zero modes motion - vibration, 
breathing - are taken into account. To get idea about the relative position of states with 
different quantum numbers - spin, isospin, fiavor, SU{3) representation, etc., it is necessary 
to calculate zero-mode quantum corrections to the energy of baryonic system. Corrections 
of this kind have been calculated first for the configurations of "hedgehog" type [§, later 
for configurations with axial symmetry ^ and for more general configurations, in SU{2) 
jnl as well as SU{3) cases [p], |3[. 

The chiral soliton approach provides the concept of nuclear matter different from that 
widely accepted as being constructed from separate nucleons, only. To find the "smoking 
gun" for this unusual concept it is necessary to find some states which cannot be made of 
separate nucleons because, e.g. of the Pauli exclusion principle. The simplest possibility is 



X 



to consider the B ^ 2 system where the Pauh principle strictly and unambiguosly forbids 
definite sets of quantum numbers for the system consisting of separate nucleons. 

Here we discuss first the SU{2) case (next Section) where supernarrow low lying 
dibaryons were predicted [jl^ and estimate their electromagnetic decay width. Further the 
SU{2>) extention of the chiral soliton model is considered and estimates for spectra of multi- 
baryons with fiavour (strangeness, charm or bottom quantum number) made previously are 
extended up to highest baryon numbers where the necessary theoretical information on mul- 
tiskyrmions is available . A simplified model for large B multiskyrmions given by rational 
maps is presented which allows us to establish the connection with the domain wall or 
bag approximation (Section 4). The technical details necessary for calculations are available 
in literature, and here some of them given in Appendices, where several statements valid for 
any chiral solitons are proved and useful expressions for SU{2) skyrmion tensors of inertia 
still lacking in the literature are presented. 



2 Narrow dibaryons below NNtt threshold 

The topological chiral solitons (skyrmions) are classical configurations of chiral fields incor- 
porated in unitary matrix U e SU{2) or SU{3) and characterized by topological, or winding 
number identified with baryon number B. The classical energy (mass) of these configura- 
tions Mci is found usually by minimization of energy functional depending on chiral fields. 
As any extended object skyrmions possess also other characteristics like moments of inertia 
8 (tensors of inertia in general case, see Appendix A) , mean square radii of mass and baryon 
number distribution, etc. The quantization of zero modes of chiral solitons allows to obtain 
the spectrum of states with different values of quantum numbers: spin, isospin, strangeness. 



etc. This approach allows for quite reasonable description of varyous properties 

of baryons, nucleons and hyperons, therefore, it is of interest to consider predictions of the 
models of this kind for baryonic systems with B > 2. The energy of SU{2) quantized states 
with axial symmetry can be presented as |10[ 

+ , J(J+1) , J| / 63 „2 03 



Here / and J are isospin and spin of the system, J3 is body fixed 3-d component of the 
angular momentum which can be considered as an additional internal quantum number of 
the system. B = n - azimuthal winding number for the lowest energy axially symmetrical 
configurations. This formula, being obtained rigorously from the model lagrangian [P, [T0[] , 
is very transparent in its physical meaning. The technical details beginning with known 
lagrangian of the Skyrme model, expressions for Md, tensors of inertia and some other 
formulas can be found in Appendix A. 

The (generalized) axial symmetry of the configuration with B ^ 2 leads to definite 
constraint on the body-fixed 3-d components of the isospin and angular momentum: J3 = 
-n/3 = -nL |10[. As a consequence of this, the states with / = 1 and J = 0, or 7 = 0, J = 1, 
and also I = j = 1 should have Ig-'^ = J^^ = L = 0. Therefore, the last term in (1) proportional 



to Jg-^^ is absent in these cases. Since the parity of configuration equals to P = (-1)^ |]T^, all 
states mentioned above have positive parity. The state with / = 0, J = 2 can have I^^ = Jg-^ = 
as well as, e.g., I^^ = L = I, J^^ = -2. At large B by special reasons (see Appendix A) also 
only the first two terms in (1), ~ /(/ + 1) and ~ J{J + 1) are important in quantum correction 
to the energy. 



As it was noted long ago the quantum correction for the deuteron- hke state 
with / = 0, J = 1, = 1/0. j{B = 2) is by about - 30Mev smaUer than that of the "quasi- 
deuteron" state with 7 = 1,7 = 0, i?™* = l/9/(i? = 2). This takes place for all known 
variants of the model, without any tuning of parameters, therefore, it can be considered as 
intrinsic property of chiral soliton models originating from effective field theories. Further 
investigations of nonzero modes of two-nucleon system, not all but many of them, have shown 



that the binding energy of the deuteron can be reduced to 6M ev [|Tg if it is considered as 
a difference between states with deuteron and quasideuteron quantum numbers. Previously 
and here we consider the differences of energies of quantized states because they are free of 
many uncertainties due, e.g., to unknown loop corrections to the masses of skyrmions (see 
TT] , |T^ and discussions below) . 



According to expression (1) dibaryons are predicted decoupled from 2- nucleon channel 
as a consequence of the Pauli principle [jl^. For example, the state with isospin / = j = i. 



positive parity and the energy below the threshold for the decay into NNtt with E'^* = 
l/Qj{B = 2) + 1/6/(5 = 2). This dibaryon cannot be seen in nucleon-nucleon interactions 
directly, but can be observed in reaction NN NNjj, where one photon is necessary to 
produce D and the second one appears from the decay of D: e.g., pp £)++7 pp^i- The 
chiral soliton models predict the state D with isospin / = j = l at the energy about 50 -60 Mew 
above the NN threshold fl^ . 

In the paper fl^ it was shown that the states with even sum / + J (O, 2, etc.) and 
positive parity are forbidden by the Finkelstein- Rubinstein type constraints which appear 
as a consequence of requirement that the configuration can be presented as a system of 
two unit hedgehogs at large relative distances, and these unit skyrmions possess fermionic 
properties. It means, that the configurations which cannot be considered as consisting of 



two nucleons, were ignored in [|Ty]. Opposite to this, we abandoned this requirement ||14|| . 
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It should be noted also that the state with / = o, J = 2 which was forbidden in 
be in fact the state of two nucleons and should not be forbidden by FR - constraint. 
Therefore, this particular case should be analyzed more carefully. 

It is possible to estimate the width of the radiative decay D NNj. Electromagnetic 
nucleon formfactors can be described quite well within Skyrme soliton model in wide interval 



of momentum transfers [^, reasonable agreement with data takes place for deuteron and 
27V - system |T^, therefore, one can expect reasonable predictions for systems with greater 
baryon numbers or with unusual properties. The dimensional estimate of narrow dibaryon 
decay width was made in providing the lower bound for the decay width of few eV. To 
make more realistic estimate one can consider magnetic type transition D NNj or dj. The 
amplitude of the direct process due to magnetic dipole transition can be written as 

Mo^NN-t = ie p.D^NN '^iklFik'^f (j)\4i2 (2) 

where jl is the value of the transition magnetic moment, we assume that /i is of the order of 
Fik = e^gfe - euqi - the electromagnetic field strength, ^f, 0i and 02 are the wave functions 
of the dibaryon and nucleons. For the width of such direct decay we obtain then 

Tu^NN, = aAM^^^^iA/My/' (3) 

which is numerically less than 0.1 for ^ fip - — 4.7/(2MAr), A = Md - 2M is the en- 
ergy release, or the maximal energy of emitted photon. This estimate agrees with that 
made previously [^|, but the final state interaction could increase it by several orders of 
magnitude. 



To take it into account roughly one should consider the transition D ^ d' where d' 
is spin zero quasideuteron, or D+ d. At this point the important statement is that the 
isovector magnetic transition operator for any skyrmion is connected simply with its mixed, 



or interference tensor of inertia Q'X- This statement, known for some particular cases P, ^ 



is proved in Appendix B for arbitrary skyrmions and for any type of chiral soliton models: 



= -ii?''^"(A)e;f Of (A'), (4) 

where i?"^ = Dlj = Tr{A^T°-AT^)/2 and Of are the final rotation matrices, a is isotopical (octet 
in SU{3)) index, and for electromagnetic interaction we should take a = 3. 9j"* is presented 
in Appendix A. 

For configurations with generalized axial symmetry, as well as for several known 
multiskyrmions, only diagonal elements of 9"* are different from zero, moreover, only the 
(33) component remains in the case of axial symmetry, and we have 

ii^ = -l^R'\A)&^:,'OUA'), (5) 

9™* = 2933 = 14.8 Gew^^ for B ^ 2 and accepted values of model parameters, see also Table 1 
below. To get numerical values of the transition magnetic moments one should calculate the 
matrix elements of rotation matrices between the wave functions of initial and final states 



which are equal in terms of final rotation matrices Dj,^ ^, see e.g. 



and we have for D state I = J = 1, L = 0, for the final d' state also / = 1, and J = 0. Since 
i?33 ^ jji^ i]^Q isotopical part of the matrix element for D ^ d' transition is proportional to 

< Dl^Dl^Dl^ j Dl^Dl^Dl^dv = Ci\o';ij3 4o;i,o/3- (7) 

The Clebsch-Gordan coefficient Ci^'oa.o = 0? therefore D d' transition magnetic moment 
equals to zero for all states including £>++ and 13°, not only for D+ d'+ which is trivial, 
and this is a consequence of symmetry property of rotator wave function with L = 0. 

For the transition D+ dj the isotopical part of the matrix element differs from zero, 
< dI qDIqDqq >= 1/3, but the angular momentum part proportional to < Dj^qDoo^j^o > equals 
to zero again. However, the decay 1?+ np is possible as a result of isospin violation in the 
second order in electromagnetic interaction, due to virtual emission and reabsorbtion of the 
photon and due to isospin violation by the mass difference of u and d quarks. The order of 
magnitude estimate of the width of this decay due to the virtual electromagnetic process is 



,2 



M 



z,_p„ i^aVvA/M (8) 
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which is about ^ iKeV. It should be noted here that for the components of D with charge +2 
or the decay into pp or nn final states is forbidden strictly, due to the rigorous conservation 
of angular momentum and the Pauli principle. 

For transition pp^, D° nn-f as well as D+ (pn)/=i7 the isoscalar magnetic 

momentum operator gives nonzero contribution. The corresponding matrix element 

Mo^d'^ = ie eMF,k^f^'''^ (9) 

The approximate relation takes place for rational map parameterization: 

A§-^3^^, (10) 



< Tq > is mean square radius of i3-number distribution. (10) coincides with result of 



for 5 = 1 and is close to result of [0 for B ^ 2. The derivation of (10) valid for rational 
maps parametrization of skyrmions, will be given elsewhere. The coefficient before J3 in (10) 
depends remarkably weak on baryon number, as can be established from Table 1. However, 
numerically (10) gives about twice smaller result for B = 1 for parameters we take here, than 
in As a result we have: 

For the decay width we obtain then 

I D^d'j — a (iZ) 

Numerically, flD~.d' ^ 0.35 Gew^^ and from (12) Tn^d'^ ^ 0.3 Kev {A/60 Mev f. The same esti- 
mate is valid for the decay rate of D+ — > npj with np-system in isospin / = 1 state. 

The experimental evidence for the existence of narrow dibaryon D in reaction pp 
pp'j^ has been obtained in Dubna [^, although these data have not been confirmed in 



the Uppsala bremsstrahlung experiment [^. Even more clear indications for the existence 
of low lying dibaryons have been obtained in experiment at Moscow meson factory in the 
reaction pd pX [|2^ . The checking and confirmation of these results in its importance is 
comparable with a discovery of new elementary particle. The absence of such states would 
provide definite restrictions on applicability of the chiral soliton approach and effective field 
theories. 

It should be noted that within the model there is a problem of the lowest state with 
I = J = which should be lower than the deuteron-like state. Therefore, deuteron should 
decay into this (0,0) state and a photon, but two-nucleon system in singlet ^6*0 state could 
not decay because 0^0 transition is forbidden for electromagnetic interaction. The loop 
corrections to the energy of states, or Casimir energy [^, are different for states which 
can go over into two nucleons, and for states which cannot. Their contribution can change 
the relative position of these states and shift the (0,0) state above the deuteron, but very 
nontrivial calculation should be made to check this. 

Some low-lying states with strangeness are also predicted, which cannot decay strongly 



due to the parity and isospin conservation in strong interactions [14]. For example, the 



dibaryon with strangeness S = -2, / = 0, J = 1 and positive parity has energy by O.lTGev 



above AA threshold ||2J], and it cannot decay into two A-hyperons because of the Pauli 
principle, and into AAtt final state by isospin conservation. Therefore, the width of electro- 
magnetic decay of such state should be not more than few tenths of Kev. It is, of course, a 
special case. Other possible states with flavour s, c or 6 will be discussed in the next section. 

The masses of neutron rich light nuclides, such as tetra-neutron, sexta-neutron, etc. 
can be estimated using formula (1). For multineutron state with / = B/2 the rotation energy 
j^rot ^ ij(i3 + 2)/(89/), and such nuclides are predicted well above the threshold for the strong 
decay into final nucleons. With increasing baryon numbers the energies of neutron rich 
states with fixed difference N - Z become lower, so, their width can be quite small. The 
mass difference of states with isospin / and ground states with / = (for even b) equals to 
AE{B,I) = /(/+l)/(2e/,B). For such pairs of nuclei as ^Li-^Be, ^^B-^^C and ^^N-^^O it equals 
to AE{B, 1) = 1/8/, B and decreases with increasing B, i.e. atomic number, both theoretically 
(see Table 1. below) and according to data. For B = 16 this difference equals to 10.9 Mew in 
comparison with theoretical value of 15.8 Mew which is not bad for such a crude model. 



3 Flavoured multibaryons 



Another characteristic prediction is that of multibaryons with different values of flavours, 
such as strangeness, charm or bottom quantum numbers. The bound state approach of 
multiskyrmions with different flavours is the adequate method to calculate the binding 
energies of states with quantum numbers s, c or b. The so called rigid oscillator model is 



the most transparent and controlable version of this method [25]. The references to pioneer 
papers can be found also in ||2^. The binding energies of flavoured states are predicted 
smaller than binding energies of ordinary nuclei - for strangeness quantum numbers, and 
greater - for charm or bottom quantum numbers. Here I present the main results on flavoured 
multibaryons following the papers and extended to higher values of baryon numbers. 

To quantize the solitons in SU{3) conflguration space, in the spirit of the bound state 
approach to the description of strangeness, we considered the collective coordinates motion 
of the meson flelds incorporated into the matrix U e SU{3), see Appendix A: 



[/(r, t) = R{t)Uo{0{t)r)R\t), R{t) = A{t)S{t), 



(13) 



where Uq is the SU{2) soliton embedded into SU{i) in the usual way (into the left upper cor- 
ner), A{t) e SU{2) describes SU{2) rotations, S{t) e S'C/(3) describes rotations in the "strange", 
"charm" or "bottom" directions, and 0{t) describes rigid rotations in real space. 



S{t)^exp(iV{t)), V{t)= Da{t)\a, 



(14) 



a=4, 



\a are Gell-Mann matrices of the (u,d, s), {u,d,c) or {u,d,b) SU{3) groups. The {u,d,c) and 
(m, d, h) SU (3) groups are quite analogous to the (m, d, s) one. For the (m, d, c) group a simple 
redeflniton of hypercharge should be made. For the {u,d,s) group, = {K+ + K^)/\/2, 
D5 = i{K+ - K-)/V2, etc., for the (w,d,c) group = + Z)")/V2, etc. 

The angular velocities of the isospin rotations are deflned in the standard way: 
A^A = -iujT/2. We shall not consider here the usual space rotations explicitly because the 
corresponding moments of inertia for baryonic systems {BS) are much greater than isospin 
moments of inertia, see Table 1., and for lowest possible values of angular momentum J the 
corresponding quantum correction is either exactly zero (for even b), or small. 

The fleld D is small in magnitude, at least, of order l/^/Nl, where Nc is the number 
of colours in QCD. Therefore, an expansion of the matrix S m. D can be made safely. To the 
lowest order in fleld D the Lagrangian of the model (^1) can be written as 



L = -M, 



cl.B 



fB(Fi 



D^D 



(15) 



m% ^ {Fl/F^)m% - ml. Here and below D is the doublet K+, K° {d", L>-, or B+, B°). ep is 
the moment of inertia for the rotation into the "flavour" direction (F = s, c, or b, the index 
c denotes the charm quantum number, except in Nc): 



Of^b = ^ / (1 - c/) 



d-'r. 



(16) 



where / is the proflle function of skyrmion, Fd is flavour decay constant, i.e. decay constant 
of kaon, or D- meson, or B-meson, 



F^ 



(1 - Cf)d\ 



(17) 



The mass term contribution to static soliton energy is connected with r due to relation 
M.t. = mJr/2. The quantity enters when flavour symmetry breaking in flavour decay 
constants is taken into account: 



(18) 



It is connected with other calculated quantities via relation: 



T = 2{M'i^/Fl-e'Qfl 

where M^^'' is second order contribution into static mass of the soliton, 8f,'' is Skyrme term 
contribution into flavour moment of inertia. The contribution proportional to in (15) 
is suppressed in comparison with the term r by the small factor ~ F|,/to|,, and is more 
important for strangeness. The term proportional to NcB in (15) arises from the Wess- 
Zumino term in the action and is responsible for the difference of the excitation energies of 
strangeness and antistrangeness (flavour and antiflavour in general case) pB |. 

Following the canonical quantization procedure the Hamiltonian of the system, in- 
cluding the terms of the order of N^, takes the form 



Hn = M, 



cl.B 



46 



F.B 



.^2 



-a 



1&Qf,b 



D^D 



NrB 



86 



'F,B 



{D^n~Ii^D). (19) 



n is the momentum canonically conjugate to variable D which describes the oscillator-type 
motion of the (w, d) SU{2) soliton in SU{i) configuration space. After the diagonalization 
which can be done explicitly [^, the normal-ordered Hamiltonian can be written as 



Hb = Mci.B + 0JF,Ba^a + LUp^Bb^'b + 0{1/N^), 



(20) 



with flt, &t being the operators of creation of strangeness, i.e., antikaons, and antistrangeness 
(flavour and antiflavour) quantum number, ujf,b and u>f,b being the frequences of flavour 



(antiflavour) excitations. D and n are connected with a and b in the following way ||25|| : 



= ^NcBKF,B{b" - a^')/(2*) 



(21) 



16(m2,rB 



{Fl~ F'^)fBQF,B)/{NcBfY/^. For the lowest states the values of D 
+ NlB^Y^'\ and increase, with increasing flavour number |F| 



with = [1 
are small: D - 

like (2|F| + 1)^/^. As was noted in |^5[, deviations of the field D from the vacuum decrease 
with increasing mass mu, as well as with increasing number of colours Nc, and the method 
works for any mo (and also for charm and bottom quantum numbers). 



ujF,B = N,B{nF,B - l)/(86f,B), u;f,b = N,B{hf..b + 1)/(86f,b). 



(22) 



As was observed in [^, the difference - ^f,b = NcB/{AQf,b) coincides, to the leading 



order in Nc with the expression obtained in the collective coordinates approach ||24 



The fiavor symmetry breaking {FSB) in the fiavour decay constants, i.e. the fact that 
Fk/F^ ~ 1.22 and Fd/F^ = 1.7 ± 0.2 (we take Fd/F^ = 1.5 and Fb/F^ = 2) leads to the increase 
of the fiavour excitation frequences, in better agreement with data for charm and bottom. 



It also leads to some increase of the binding energies of BS ^6 . 
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Table 1. Characteristics of the bound states of skyrmions with baryon numbers up to i? = 22. The classical 

mass of solitons Md is in GeV , moments of inertia Of, B/, 0/,3, and Oj as well as < ro >, F and F - in GeV~^, 
the excitation frequencies for flavour F, cop in GeV. < ro >= with Oj defines the value of isoscalar 

magnetic momentum of multiskyrmion. For larger baryon numbers, beginning with B = 9, calculations are 
made using rational maps (RM) ansatz. For B = 32 it was assumed that the ratio I/B"^ = 1.28 as for 
RM B = 22 skyrmion. The external parameters of the model are Fj^ = 186 MeV, e = 4.12. The accuracy of 
calculations is better than 1% for the masses and few % for other quantities. 



6 J shown in Table 1 is 1/3 of the trace of corresponding tensor of inertia, see Appendix 
A. As it can be seen from Table 1 the flavour excitation energies increase again for B = 22, 
and the important property of binding becomes weaker for largest B. It can be, however, the 
artefact of the RM approximation discussed in the next Section. In particular, for rational 
maps solitons with i? > 9 we take as moment of inertia 6/ and 9/,3 1/3 of the trace of 
corresponding tensor of inertia, see Appendix A. 

For large value Fd/Ft, = pd and mass mo, the following approximate formula for the 
flavor excitation frequences can be obtained: 

. / Qf%\ N,B 

with m% = m\,+F^VB/TB- It is clear from (23) that, first, w's are smaller than masses of mesons 
mc, i.e. the binding takes place always, and it is to large degree due to the contribution 
of the Skyrme term into the flavour inertia Of.''. When po ^ oo, u>f ^ mo- Since the ratios 
Tb/Tb decreases with increasing B and &f,b/^b increases when B increases from 1 to 4 - 7, 
the energies u>f,b decrease for these B-numbers, therefore, it leads to increase of binding of 
flavoured mesons by SU{2) solitons with increasing B up to S 4 - 7. However, for B = 22 
and 32 the ratio Qf.b/^b is smaller than for i? = 1, and, indeed, the w's are the same and 
even larger than for B = 1. 
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Table 2. The binding energy differences Aes^c,b are the changes of binding energies of lowest BS with 
flavour s, c or 6 and isospin I = Tr + \F\/2 in comparison with the usual u, d nuclei, for the flavour numbers 
S = -1, -2, c = 1, 2, 6 = -1 and -2 (see Eq. (24)). The SU{3) multiplets are (p, q) = (0, 3B/2) for even B 
and {p, q) = (1, {3B - l)/2) for odd B. 



The binding energies differences between flavoured multibaryons and ordinary nuclei 
in the rigid osciUator approximation are given by the formula: 



AeB,F = \F\ 



3(/tf,i - 1) TrinpM - 1) {\F\ + 2){kf,b - 1) 
dtHi^p^QF,! '^kf^bQf.b 



(24) 



and the lowest SU{3) multiplets are considered with isospin of flavourless component Tr ^ 
for even B and Tr = 1/2 for odd B. This formula is correct for |F| = 1 and for any \F\ if the 
baryon number is large enough to ensure the isospin balance. 

The values of Ae shown in Table 2. should be considered as an estimate. They 
illustrate the restricted possibilities of RM approximation for large B multiskyrmions. 

Isosinglet BS, in particular those with \F\ = B are of special interest. As it was 
argued in such states do not belong to the lowest possible SU (3) irreps, they should have 
Tr = \F\/2. It makes sense to calculate the difference of binding energy of such state and 
the minimal state (^'""^g'"™ with zero flavour which we identify with usual nucleus (ground 
state): 



AeB,F = |i^| 



UJF,1 — I^F,B 
1 



3{kf,i 



1) _^{jF\+2)i>,F,B-l) 



[\F\{\F\ + 2)/4 - r;"™(r™" + 1)] 



where T™*" = 0, or 1/2 as before. 



(25) 
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Table 3. The binding energies differences of lowest flavoured BS with isospin / = and the ground state 
with the same value of B and 7 = or 7=1/2. The first 3 columns are for |F| = 1, the next 3 columns - for 
\F\ = 2, and the next 3 - for \F\ = 3. The state with the value of flavour \F\ belongs to the SU{3) multiplet 
with Tr — \F\/2. In the last column the binding energies differences are shown for the isoscalar electrically 
neutral states with S — —B. For |F| > 3 all estimates are very approximate. 

According to Table 3 the total binding energy of state e.g. with i? = 22 and S = -2 
is by 73 Mew smaUer than that of nucleus A = 22, so it should be well bound. The model 
used here is too crude for large values of flavour, and results obtained can be used only 
for illustration and as a starting point for further investigations. Results similar to those 



described in this section are obtained also in other versions of the model |27], in particular 



in the quark-meson soliton model [28|. For baryon numbers 73 = 3,4 estimates of spectra of 
baryonic systems with charm quantum number were made in |2^] within conventional quark 
model. They are in fair agreement with ours. 

In the channel with 73 = 2 the near-threshold state with strangeness S = -1 was 
observed long ago in the reaction pAK+ [RD| and conflrmed recently in COSY experiment 



Similar near-threshold AA state was observed by KEK PS E224 collaboration [^. The 
Skyrme model explains these near-threshold states with 7? = 2, and predicts similar states 
for greater values of the 73- number. For some values of 73, beginning with 73 > 5, 6 such states 
with several units of strangeness can be stable relative to strong interactions. Due to well 
known relation between charge, isospin and hypercharge of hadrons, Q = 73 + (73 + S)/2, the 
735" with several units of strangeness can appear as negatively charged nuclear fragments. 
For even B and minimal multiplets {p,q) = (0,373/2), strangeness S = -27, and condition when 
Q = -1 fragment appears first is -1 = 5 + 73/2, or -5 = 73/2 + 1. For 73 = 6 it is S* = -4, for 
73 = 8, S = -5, etc. For odd 73 the Q = -1 state should have strangeness -5 = (73 - l)/2 + 1, 
i.e. -3, -4 and -5 for 73 = 5, 7 and 9, etc. 

The negatively charged long living nuclear fragment with mass about 7.4 Gev observed 
in NA52 CERN experiment in Pb + Pb collision at the energy of 158 A Gev can be, within 
the chiral soliton models, a fragment with 73 = 7 or 6 and strangeness S" = -4 or -5, -6. 
Confirmation of this result as well as searches for other negatively charged fragments would 
be of great importance. For charm or bottom quantum numbers the binding energies are 
greater, but to observe such states one needs considerably higher incident energies. 



4 Large B multiskyrmions from rational maps in the domain wall 
approximation 

The rational map ansatz for skyrmions proposed in and widely used now, also in present 



paper, simplified considerably the treatment of multiskyrmions, and, at the same time, it 
leads to the picture of multibaryon system at large 73 which is, probably, incompatible with a 
picture for ordinary nuclei. To clarify this point, we consider here large 73 multiskyrmions in 
some kind of a toy model - in a domain wall approximation, which gives, however, quite good 
numerical results for known RM multiskyrmions, except i3 = 1,2. The energy of skyrmion 



for rational map ansatz [15| in universal units 37r^F^/e is: 



1 



M^^ll Aj,r^f'^ + 2734(/'2 + 1) +x4 \dr (26) 



in \ r 



The coefficient An = 2{N-1)/N for symmetry group SU{N) [Q. Tfie quantity I for 5L/(2)-case 
is given in Appendix A, tlie inequaiity talces place I > B^. Direct numericai caicuiations liave 
sliown, and tfie anaiytical treatment fiere supports, tfiat at iarge B and, lience, iarge I muiti- 
slcyrmion iooks like a spfiericai ball witti profile equal to F = tt inside and F = outside. The 
energy and B-number density of this configuration is concentrated at its boundary, similar 



to the domain walls system considered in in connection with cosmological problems. 

Consider such large B skyrmion within the "inclined step" approximation. Let W 
be the width of the step, and tq - the radius of the skyrmion where the profile / = n/2. 
/ = 7r/2 - (r - ro)n/W for - W/2 < r < Vo + W/2. Note that this approximation describes the 



usual domain wall energy |^ with accuracy 9%. 

We write the energy in terms of W, , then minimize it with respect to both of these 
parameters, and find the minimal value of energy. 



M{W,ro) 



This gives 



W„ 



1/2 



and, after minimization, r?.. 



B + ZI/{Srl)_ 
In dimensional units 



(27) 



(28) 



VW(8^^v)- In dimensional units ro = (6X/AA,)i/''/(F^e). Since 
I > B^, the radius of minimized configuration grows as Vb, at least. Wmm = tt, i.e. it does 
not depend on B for any SU{N). The energy 



{2B + v/3AArX/2)/3 



(29) 



For SU{2) model An = 1 and the energy Mrmn = (25 + y3x72)/3 should be compared with the 
lower bound Mlb = (2B + \/7)/3. The formula gives the numbers for 5 = 3, 22 in remarkably 
good agreement with calculation within RM approximation, within 2-3% 0. 

It is not difficult to calculate the corrections to these expressions, of relative order 
1/B, l/B\...: 



/3 = ^V(12S), 7^ 



M(W^,ro)~ 
(27r2 + 17)/\/24X. 



W 



{B + A^rD + W 



B{\ + P) 



[2B(l+/3/2) + v/3l78(l+7/2)] 



(30) 



(31) 



However, the first order correction in W does not improve the description of masses, and 
summation of all terms seems to be required. 

So, we see that a very simple approximation provides a confirmation of a picture from 
numerical calculation of RM skyrmions as a two-phase object, a spherical ball with profile 
/ = TT inside and / = outside, and a fixed width envelope with fixed surface energy density, 
pj^i = M/(47rr§) ~ {2B + ^ZI / 2) / (\2t: rl). At large B pM Const, but average mass density over 
the volume 0. 

Consider also the influence of the mass term which gives the contribution 



M.t. =fn J r-^(l - cosF)dr, 



(32) 



in = 8m^/(37rF^e^). For strangeness, charm, or bottom the masses m^, mo or should be 
inserted instead of m^r. In the "inclined step" approximation we obtain: 



M.t. 



(33) 



In view of this structure of the mass term it makes no influence on the width of the step W 
in lowest order, but the dimension of the soliton Tq becomes smaller: 



ro^To-m — . (34) 

As it was expected from general grounds, dimensions of the soliton decrease with increasing 
TO. However, even for large value of m the structure of multiskyrmion at large B remains 
the same: the chiral symmetry broken phase inside of the spherical wall where the main 
contribution to the mass and topological charge is concentrated. The behaviour of the 
energy density for S = 22 at different values of ^ is shown in Fig. The value of the mass 
density inside of the ball is defined completely by the mass term with 1 - c/ = 2. The baryon 
number density distribution is quite similar, with only difference that inside the bag it equals 
to zero. It follows from these results that RM approximated multiskyrmions cannot model 
real nuclei at large B, probably B > 12-20, and configurations like skyrmion crystals may 
be more valid for this purpose. 

Besides the simple one-shell configurations considered in [^, |15| and here, multi- 
shell configurations can be of interest. Some examples of two-shell configurations with 
B = 12, 13, 14 have been considered recently [^^. The profile / = 27r at r = for such configu- 



rations and decreases to / = for r oo. We can also model such two-shell configuration in 
the domain-wall, or spherical bag approximation with a result 

M~ (2Bi + v/3X72)/3+(2B2 + V3V2)/3, (35) 

with total baryon number B = Bi + B2. The profile / decreases from 2tt to n in the first 
shell, and from tt to in the second. The radii of both shells should satisfy the condition 
''o^^ > '^0^'' + so external shell should be large enough, with baryon number B2 of several 
tens, at least. Since the ratio I/B'^ is greater for smaller B, the energy (35) is greater than the 



energy of one-shell configuration considered before. Calculations performed in ||3^ also did 
not give results better than for one-shell configurations. However, more refined consideration 
would be of interest. Observation concerning the structure of large B multiskyrmions made 
here can be useful in view of possible cosmological applications of Skyrme-type models. 



5 Concluding remarks 

Here we have restricted ourselves with the Skyrme model and its straightforward extentions. 
However, many of the result are valid in other variants of the model: in the model with 
solitons stabilized by the explicit vector (w) meson, or stabilized by the baryon number 
density squared, in the chiral perturbation theory, etc, see discussion in [jl^b. The B = 2 



torus-like configuration has been obtained within these models, as well as in the chiral quark- 
meson model , and it would be of interest to check if there are also multiskyrmions with 
B>3. 

We did not discuss a special classes of SU{3) skyrmions, S0{3) solitons and the prob- 
lems of their observation, as well as SU{3) skyrmion molecules. The discussion of these 
topics can be found in several papers of IT^, O. Some new solutions which are not SU{2) 



embeddings in SU{3) or SU{n) have been found in . 

To conclude, the study of some processes, also at intermediate energies which, to some 
extent, are out of fashion now, can provide a very important check of fundamental principles 
and concepts of the elementary particles theory including the confinement of quarks and 
gluons. The confirmation of the predictions of the chiral soliton approach would provide 



qualitatively new understanding of the origin of nuclear forces. If the existence of low energy 
radiatively decaying dibaryons is reliably established, it will change the long standing beliefe 
that nuclear matter fragments should consist necessarily of separate nucleons bound by their 
interactions. Therefore, the confirmation and checking of the results of experiments on 
dibaryons production, as well as of production of fragments of flavoured matter is extremely 
important. It would be possible at accelerators of moderate energies, like COSY (Juelich, 
FRG), KEK (Japan), Moscow meson factory (Troitsk, Russia), ITEP (Moscow), and some 
others. The production of multistrange states, as well as states with charm or bottom 
quantum numbers, is possible in heavy ion collisions, and also on accelerators like Japan 
Hadron Facility to be built in the near future. 

The multiple flavour production realized in the production of flavored multibaryons 
possible, e.g., in heavy ion collisions, demands higher energy, of course, but multiple in- 
teraction processes and normal Fermi motion of nucleons inside of nuclei make effective 



thresholds much lower ||3^. Studies of such flavoured multibaryons production would allow 
more complete and reliable checking of the model predictions. 

We note finally that the low energy dibaryons have been obtained recently in 
using the quantization procedure different from our. 

The work is supported by RFBR grant 01-02-16615, UK PPARC grant PPA/V/ 
S/1999/00004 and presented in part at the International seminar Quarks-2000, Pushkin, 
Russia, May 2000. 

Appendix A. Inertia tensors of multiskyrmions. 

The lagrangian density of the SU{2) Skyrme model is given by 

p2 1 F2 2 

^ = -Tt^"^^^- + 32^^^<^ + ^^KC/ + U^- 2), (.41) 

Lfj, = df_,UW is left chiral derivative, = zL^^feTfe, Tk are Pauli matrices. G^i. = d^L„ - d^L^^ - 
strengh of the chiral field. The Wess-Zumino term present in the action has been discussed 



in details in [|T^, and we shall omit this discussion here. 

First we give the expression for the energy of SU{2) - skyrmion as a function of profile 
F and unit vector fi, which is especially useful in some cases. Using definition U ^ Cf + isfnf, 
and relation 

Lf,,kL^,k = d^fdvf + sjdf.fid^n, {A2) 



we obtain 



+ PAi.t.>d^r. (AS) 



For the Ansatz based on rational maps the profile F depends only on variable r, and 
components of vector n - on angular variables 9, (f>. = {2ReR)/{l + |i?P), Uy = (2ImR)/(l + 
|i?P), = (l-|i?p)/(l+|i?p), where i? is a rational function of variable z = t5(6'/2)ea;p(i(/)) defining 
the map from S*^ S^. In this case the gradients of functions F and n are orthogonal 
(recall that dr = firdr + ngde/r + n^d^/{rse) Ur = f/r = (sec^, ses^, cg), He = {-cqc^, -ces^, se), 
= (s0, -C0, 0)) and [dfdni]'^ = /'^(9ni)^, etc. Taking into account relations 



ni[dn2dn3Y = ni[dnidn2Y, nUdnidnaY = ni[dnidn2Y , (A4) 
one can present {A3) as 



f\dn,f + s}[dn,dn2f/nl 



+ PM.t.>d^r. (A5) 



Usually the notation is introduced 



(l + |z|2) \dR\\ 2idzdz 



An J \{i + m \dz\ ; (i + |2p)2' 



{A6) 



and using the equation 



J m\ 



dn = 2 



1^ 
J (1 



2idRdR 



obtain finally 



stat 



47r J I f (/'V^ + 2s}Af) + ^ [2/'^ + sp/r 



+ r pM.t. >dr. 



{A8) 



To find the minimal energy configuration at fixed M = B one minimizes I, and then finds 
the profile F{r) by minimizing energy (AS). 

To quantize zero modes one uses ansatz U{t,r) = A{t)U{Oik{t)rk)A^t), and evident 
relation 

dtU = U = AU{r{t))A'^ + AU{r{t))A'^ + fi{t)AdiU{r{t))A\ {A^) 

where ri{t) = Oik{t)rk is body- fixed coordinate. 

Angular velocities of spatial (or orbital) rotations are introduced according to: 

ri = Oikr'k = OikOllnit) = -eum^rnn{t) 

and integration is performed in coordinate system bound to soliton (body-fixed). 

The rotation, or zero-mode energy of SU{2) skyrmions as a function of angular veloc- 
ities is 



The isotopical tensor of inertia for arbitrary SU{2) skyrmion is: 



4 e2 

For the RM ansatz trace of this tensor of inertia is 



+ ^rdiUadmb }d^r. 



(AlO) 



{All) 



{A12) 



The orbital inertia tensor gives contribution to the energy Q-^^Qa^bf^^, and for arbitrary 
configuration using the same notations is given by: 



difdkfidnif + [dffdindkri- 



-difdifdindk-n - dkfdifdindifi + s'j [{dn{f' difidkn - {dindin){dkndin)] 
This expression can be simplified for RM ansatz: 



>eiaa^kl3brarf3(fr. 



(A13) 



^ab 



{dni)^{r^Sab - Tan) - daudbfir^ 



f/ 
«2 



{dindkn){dindkfi){r^6ab - ran) - r'^{dandin){dbndin) 



It allows to obtain easily the trace of the inertia tensor 9 



J 



{AU) 



{A15) 



It is easy to establish the inequality for traces of isotopical and orbital tensors of 



inertia: 



ei - seL = ^{x-B^)J sjdr > 0, 



{AW) 



since I > B'^, Af = B. The interference (mixed) tensor of inertia which defines also the 
isovector part of the magnetic transition operator equals to: 



0int 
ab 



^ + ^Aid^fr + 4id.n)'] 



diTir 



— ^{difduf + s'fdind„n)d„ni ynkekia^iabTad^r. 



{All) 



The components of spatial angular velocities interfere with components ui, u>2, of angular 
velocities of rotation in configuration space, only. 

Numerically the components of the mixed tensor of inertia are much smaller than 
those of isotopical, or orbital tensor of inertia, except special cases of "hedgehogs" when 
_Qint = 9/ = and axially symmetrical configurations when for the 3-d components of 
inertia relations hold -©la* = nS^g = ©33 /n. 

Note, that most general formulas for tensors of inertia are presented here for the first 
time. For the case of RM configurations they differ in some details from those given in the 
literature. 



Appendix B. Electromagnetic transition operators. 

Here we prove, for completeness, in general form some statements concerning isovec- 
tor (octet in SU{3) case) vector charge and isovector magnetic momentum operator. 

There is the following connection between isovector current and isospin generator 

Vo,a = ^Tr{A^\aAXt)ll^' = Rab{A)l'J , {Bl) 
where the isospin generator in body-fixed (connected with soliton) coordinate system is 

ll^ = dU^^w, n)/dujb- {B2) 
a, 6= 1,2,3 for S'[/(2)-model, and a,h = 1,...8 for S'[/(3)-model. To prove this consider ansatz 

U = e-'"»^»/2^(t)C/oAt(t)e*«"^"/2 (B3) 

The vector Noether current is a coefficient before derivative of the probe function, d^a. In 
the lowest order in a we obtain for the chiral derivative: 

U^doU = A[U^A^{A - iaA/2)Uo -A'' {A- iaA)]A'' {BA) 

Using the definition of angular velocities of rotation in configuration space Ua, we obtain 

AU - iA^aA/2 = -'^\b{uJb + Rab{A)aa), {B5) 



where real orthogonal matrix 

Rab{A) = ^TriA^XaAXb). (B6) 

Since the dependence on a reduces to simple addition to angular velocity according to (-B5), 
formula {Bl) follows immediately. 

According to the well known relation, 

Q = B + Is/2 = B + Vo,3/2 {B7) 

the baryonic (topological) charge and 3-d component of the isospin generator contribute 
to the charge of the quantized skyrmion. 

We prove also that there is simple connection between isovector (octet for SU{3) 
model) magnetic momentum operator of the skyrmion and mixed (interference) tensor of 
inertia. Note first that the lagrangian of arbitrary chiral model, not only Skyrme model, 
because of Lorentz invariance can be presented as a sum, with some coefficients, of contri- 
butions of the type: 

£m,jv = Tr{U'^UMU'^UN - U^'dkUMU^'dkUN) , {B8) 

where M and N are some matrices. E.g., for the second order term M = N = 1. The 
contribution into rotational energy, proportional to fl, uj, which comes from the first term 
in {B8) and defines mixed or interference tensor of inertia is (see {A9) above): 

eifwaflb = J Tr{U^A^AUo - A^A)MUldkUoNfkd^r +(^M^N^, (B9) 

M = A'^MA^N = A^NA. Or, 

©If = -\^bjk j r,{t)Tr{Ul\aUQ - KWUldkU^Nd^ + (^M ^ (BIO) 

rj{t) and dk are body- fixed here. From the second term in expression (B8) we obtain for the 
spatial components of the vector current: 

= ^Tr (C/oUtA„^f/o - A^XaA)MUldkUQN + (^M ^ {Bll) 

Taking into account that A'^XaA = Rab{A)Xb., Rab = \TrA'^XaAXb and dk = Oikd\^ ^ 

= '-RabOikTr{ulXbUQ - Xb)MUldiUoN + (^M ^ {B12) 

By definition 

= \ujk j rjV.-d'r, (S13) 

/z» = ^eijkRab{A)OgkOpj J rp{t)Tr{UlXbUo - Xb)MUldqUoN + (^M ^ N^. (B14) 
Taking into account that 

we obtain the desired relation between components of the magnetic momentum operator 
and mixed tensor of inertia in the body-fixed coordinate system: 

< = -lRab{A)eirOu. {B15) 



or 



In some particular cases this relation was used previously [§, [T^ . 

For the transition matrix elements calculations it is necessary to average this expres- 
sion over wave functions of some initial and final states, see Section 2. 
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Figure caption. 

The mass density distribution of the rational map multiskyrmion with B = 22 as a 
function of the distance from center of skyrmion for different values of mass in the chiral 
symmetry breaking term . 

a) pion mass in the mass term, b) kaon mass, c) D-meson mass, the mass density is 
devided by 10. 



